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Abstract 

By considering the master equation of the totally asymmetric exclusion 
process on a one-dimensional lattice and using two types of boundary con- 
ditions (i.e. interactions), two new families of the multi-species reaction- 
diffusion processes, with particle-dependent hopping rates, are investi- 
gated. In these models (i.e. reaction-diffusion and drop-push systems), 
we have the case of distinct particles where each particle A a has its own 
intrinsic hopping rate v a . They also contain the parameters that control 
the annihilation-diffusion rates (including pair-annihilation and coagula- 
tion to the right and left). We obtain two distinct new models. It is shown 
that these models are exactly solvable in the sense of the Bethe anstaz. 
The two-particle conditional probabilities and the large-time behavior of 
such systems are also calculated. 

Keywords: driven diffusive systems (theory), quantum integrability (Bethe 
ansatz), stochastic processes (theory) 

1 Introduction 

Our understanding of non-equilibrium statistical physics is much behind than 
that of equilibrium theory. Asymmetric exclusion processes (ASEP) are paradig- 
matic models for systems which are far from equilibrium. Despite their greatly 
reduced complexity they capture various fields of physics like the study of shocks 
[1, 2], the noisy Burgers equation [3], polymers in random media, dynami- 
cal models of interface growth [4], the traffic models [5], and the kinetics of 
biopolymerization [6]. 

On the one-dimensional lattice, the simplest system is the totally asymmet- 
ric exclusion process (TASEP). In this model, each lattice site is occupied by at 
most one particle and all particles can only hop with equal rate to their right- 
neighboring site, provided this site is empty; otherwise the attempted move 
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is rejected. In [7], TASEP has been solved by introducing a master equation 
which describes the evolution equation of the particles when they are not in 
neighboring sites, and a so-called boundary condition, which specifies the situ- 
ation in which the probabilities go outside the physical regions. This happens 
when some of the particles are in adjacent sites and the non-physical probability 
terms appear in the master equation. The coordinate Bethe ansatz (BA) has 
been used to obtain the TV-particle conditional probabilities of TASEP. 

The interesting point is that one can find the new solvable models (including 
interactions besides diffusion to right- neighboring sites), with other boundary 
conditions and the same master equation. In [8], a similar technique has been 
used to solve the so-called drop-push model. In this model the particle hops 
to the next-right site, even if it is occupied. It can hop by pushing all the 
neighboring particles to their next-right sites, with a rate depending on the 
number of these particles. Some other generalizations of TASEP can be found 
in [9, 10, 11]. 

The other interesting point is the study of multi-species systems in which 
several kinds of particles move and interact (i.e. reaction-diffusion) on a lattice. 
In [12], it has been shown that more-than-one species systems are solvable in 
the sense that the S'-matrix corresponding to them is factorizable into two- 
particle S-matrices. It has been found there that the criterion for this is that 
the interactions must be such that the S'-matrix satisfies a kind of spectral 
Yang-Baxter (SYB) equation. The multi-species generalization of the reactions 
considered in [12], has been studied in [13] and the authors remarked that the 
SYB reduces to a non-spectral matrix equation. The drop-push reaction of [8], 
has been generalized to p-species in [14]. The most general totally asymmetric 
reaction-diffusion process, including the extended drop-push interactions has 
been studied in [15]. 

The other family of diffusion system is the processes in which particles can 
hop to the next-right and the next-left neighboring sites. In [7], the single-species 
model with only simple diffusion to the right and left, i.e. partially asymmet- 
ric exclusion process (PASEP), has been studied and in [16], the one-species 
PASEP drop-push model has been described. The multi-species generalizations 
of PASEP models in which each particle hops to the right by rate Dr and hops 
to the left by rate Dl, have been studied in [17]. The Dl = case of the models 
studied in [17] leads to the previous mentioned TASEP models. 

In all of the above studies, particles hop randomly in continous time on the 
integer lattice Z. A hopping event occurs independently for each particle with 
identical rate. The model of one-dimensional PASEP in which the rates of jump 
are chosen randomly at time zero and fixed for the rest of the evolution, was 
first introduced by Benjamini et al. [18] who proved the existence of a critical 
density p* depending only on the distribution of hopping rates. They also proved 
for this case the hydrodynamic limit of an associated zero-rang process (ZRP) 
which is obtained from TASEP by identifying particles with the sites of a new 
one dimensional lattice and the interparticle distance (number of empty sites 
between particles i, i + 1) as occupation number at site i of that lattice. 

This model has received renewed attention because of the occurrence of a 
condensation transition analogous to Bose-Einstein condensation and because 
of its close relationship with exclusion processes. An analogous transition also 
occurs if the rates are identical for each particle, but dependent on the lattice 
distance to the next particle, see [19] for a recent review (in terms of the ZRP) 
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and [20, 21, 22, 23, 24, 25] for the current developments. 

The generalization of TASEP (only for one-species) with constant hopping 
rates to the particle-dependent hopping rates (PDHR) has been studied in [26] . 
The authors considered the case that a particle in site Xi has a hopping rate 
Vi and used the BA to obtain in a determinant form the exact solution of 
the master equation for the conditional probabilities of TASEP with PDHR. 
They also derived a determinant expression for the time-integrated current for 
a step-function initial state. In the sense of the BA, the two-species TASEP in 
which different particles hop with different rates and fast particles stochastically 
overtake slow ones, has been considered by Karimipour in [27]. 

In this paper we are going to study the multi-species TASEP with PDHR. 
We obtain two new families of multi-species diffusion processes with PDHR. 
The first family (i.e. reaction-diffusion models) has the following reactions 

A Q — > %A a with rate v a , 
A a A fj -> A~ ( A S with rate c°f . (1) 

and the second family (i.e. drop- push models) , describes the following processes 

A Q — > %A a with rate v a , 
A a A$ ->• %A 1 A S with rate 

(2) 

where the dots indicate the other drop-push reactions with n-adjacent particles. 
In these models, we have the case of distinct particles where each particle A a 
diffuses to the right with its own intrinsic hopping rate v a . We also consider the 
annihilation-diffusion (including pair-annihilation and coagulation to the right 
and left) extension of these above models. The annihilation processes are 

A a A/3 — > %Ap with rate 8 a p, 
A a Ap — > A a % with rate "f a /3, 

A a Af) -> with rate r) a p. (3) 

We show that the reaction rates of these models (now including the annihilation 
processes) must satisfy some specific constraints and these models are exactly 
solvable in the sense of the BA provided a non-spectral matrix equation is 
satisfied. 

The paper is organized as follows. We first introduce two types of boundary 
conditions that can be used to obtain two new families of p-species reaction- 
diffusion processes with PDHR in section 2. We use the type 1 boundary condi- 
tion and the type 2 boundary condition, in terms of two p 2 x p 2 matrices c and 
b respectively, to generalize p-species reaction-diffusion and drop-push systems 
with identical hopping rates into PDHR systems. In section 3 we consider the 
annihilation processes extension of these new models. In section 4 we investi- 
gate the BA solution of the new models and discuss under what conditions, one 
can use the BA to find exact solutions. We show that the matrix c(b) (a ver- 
sion of c(b) that constructs the S'-matrix and determines the coefficient of BA) 
must satisfy a non-spectral matrix equation. Then in section 5 we show that 
for p = 2 (for both models), the specific class of parameters, which corresponds 
to the reactions rates together with the annihilation-diffusion rates, satisfies a 
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non-spectral matrix equation. Finally in section 6 we study the conditional 
probabilities of these models and specially for two-particle systems of section 5, 
we obtain the exact expressions and the large-time behavior of such systems. 

2 TASEP-Generalization 

2.1 Boundary condition 

Consider ap-species system with particles A\, A2, ■ ■ ■ , A p . The basic objects we 
are interested in are the probabilities P ai -a N (xi,--- ,XN',t) for finding at time 
t the particle of type a\ at site xi, particle of type 012 at site X2, etc. We take the 
physical region of coordinates as x\ < X2 < ... < xn- The master equation for 
a totally asymmetric exclusion process, with particle-dependent hopping rates, 
is 

N 

^ ^ ^CLi Pa-L---a.N {x\ j ' ' " ) Xi—i, X{ 1, Xi-\.\ , • • • , Xn, t) 
i=l 

N 

(x lt ■■■ ,x N ;t) (4) 

This equation describes a collection of TV-particle that the ai-th particle drifts 
to the next-right site by rate v ai where v ai s are finite non-zero real numbers. 
This master equation is only valid for Xi < x i+i — 1. For Xi = x i+ i — 1, there will 
be some terms with Xi = Xi+\ in the right-hand side of (4), which are out of the 
physical region. But one can assume that (4) is valid for all the physical regions 
Xi < x i+ i by imposing certain boundary conditions for Xi — x i+i . Different 
boundary conditions introduce different interactions for particles. Following the 
same argument given in [17], it can be easily seen that the master equation (4) 
leads to the following relation for two-particle probabilities 

d 

^ ' ^ ' Puict2 i x l i x 2 j t) ^ ] v ct2 Paia2 (^) ^) ~ ^ ' v ai Pg\ai {x, X -\- 1; t) 

X 2 X±<X 2 X X 

(5) 

This equation leads us to take P ai a 2 { x i x j t) as the following two types of bound- 
ary conditions 

v a2 P aia2 {x,x) =^ci\ a 2 2 Pp 1 /3 2 (x,x + 1) type 1 (6) 



v a2 P aiCt2 (x,x) = b'i^PfJif^ - ^ 2 ( 7 ) 



stands for ((iifa) and b and c are p 2 x p 2 matrices which determine the 
interactions. In the probabilities appear in (6) and (7), we have suppressed all 
the other coordinates and the time t for simplicity. We consider the processes 
in which the number of particles is constant in time, in other words we exclude 
the creation and annihilation processes (in fact in this step) . Therefor if we sum 
(5) over a\ and c*2, the left-hand side becomes zero and results in 



dt 



,(xi, ■ ■ ■ ,x N ;t) 
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EE v ai p aia2 (x,x + 1) + J2 E E c tt o*. x+i) = o. (8) 

x a x 3 \ a / 

in which (6) has been used. Clearly (8) gives 

E c "i"2 = v f )l constraint for type 1. (9) 

In the same way, if we use (7), the sum over the elements of each column of 
matrix b results in 

E ^aiS = v 0i constraint for type 2. (10) 

a 

Remember that the constraints (9) and (10) are the consequence of the conser- 
vation of probabilities so they can be interpreted as the probability conservation 
equations. 

2.2 P-species reaction-diffusion systems with PDHR 

Following the same steps as [17], we first consider P aiCC2 (x,x + 1). Using (4) 
and the type 1 boundary condition (6), it is found 

P aiCt2 (x,x + l) = v ai P aia2 (x - l,x + 1) + ^ cP 1 f a 2 2 Pf !l fi 2 (x,x + 1) 

0^a 

— Va 2 P ai a 2 (x, X + 1) — ^""^ c glg^ Pg 1 a 2 ( x j X + 1) . (11) 

0^a 

in which we have used (9). It can be written as (in fact the diagonal elements 
are not reaction rates) 

„aia 2 _ _ „ai<*2 (10) 

It is seen that the evolution equation (11) describes the following two-particle 
interactions 

A a — > %A a with rate v a , 
A a A g -> A y A s with rate c"f . (13) 

It is simple to show that our formalism is consistent for morc-than-two particle 
systems. In this model, in addition to interactions, we have the case of distinct 
particles where each particle A a has its own intrinsic hopping rate v a . It is 
more general than TASEP model has been studied in [13], in which all particles 
have the equal rate. 

2.3 P-species drop-push systems with PDHR 

In the same way, we first consider P aiCt2 (x,x + 1) by using (4) and the type 2 
boundary condition (7). The result is 

P aia2 (x,x + 1) = v ai P aiCt2 {x-l,x + l) + ^b a \^ 2 P hB2 {x-l,x) 
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- v aa P aiaa (x,x + l)-Y,Vfc%P a i°,fr* + l)- ( 14 ) 

a 

in the above equation we have directly used constraint (10) as 

£« 2 =*V- (15) 



so in agreement with (14), the allowed processes are 

A a % — > ®A a with rate v a , 
A a Ap9 -> ®A 7 As with rate b"$ . (16) 

In this model any particle can hop to the right site, with rate depends on the type 
of the particle, if that site is empty. If the right site is occupied, the left particle 
can still hop to that site by pushing the right one, but in the mean time there 
is a probability that the types of the particles change. It is important to note 
that all elements of b (including the diagonal elements) should be nonnegative 
(since they are rates). We will consider more-than-two particle systems when 
we introduce the type 2 model. 



3 Annihilation-diffusion processes 

The annihilation process is a diffusion-limited reaction-diffusion process. In 
this process, the particles annihilate pairwise or coagulate to the right and left 
whenever they meet each other. Now we add the annihilation-diffusion to the 
previous reactions of both models. Note that the annihilations appear only in 
the sink terms of the evolution equation, as if we consider the initial state with 
n particles, no annihilation processes can lead to a n-particle state at any later 
time. So we do not have the conservation of probabilities (particles) and one 
can enter a sink term, A^ 2 , into the conservation equations (constraints) as 

(17) 

a 

Y, b tt=^~^ (18) 

a 

and if we use these modified constraints in the calculation of P ai a 2 {x,x + l), we 
find the same equation as (11) or (14) except extra term X aia2 P aia2 (x, x + 1) 
which is added to the sink terms. So A ai a 2 is the sum of the rates of all 
annihilation processes with initial state (ai,^) and therefor it is a positive 
quantity. These processes are 

A a A/3 — > $Ap with rate 5 a p, 

A a Ap — > A a $ with rate jap, 

A a Afs -> with rate r) a p. (19) 
and the relation between the above rates is 

K[3 = 5 afj + "fa/3 + Va0- (20) 
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One should note that the p-species reaction-diffusion and drop-push systems, 
including the annihilations, have three kinds of processes. The first one is the 
pure diffusion with PDHR which occurs when a particle is adjacent to a hole, 
the second and third ones are the reactions and the annihilations respectively, 
which occur when two particles are adjacent to each other. The master equation 
for the diffusion with PDHR is (4) and for this case the master equation does 
not contain anything about both the reactions and the annihilations. The effect 
of reactions has been coded in the boundary condition which is in fact the part 
of the master equation. Also the effect of annihilations has been coded in the 
modified constraint (17) or (18). 

4 The Bethe ansatz solution 
4.1 Solvability criteria 

To solve the master equation (4) with two types of boundary conditions, we 
consider the modified Bethe ansatz of the form [27] 

N 

P ail ..., aw (x;i) = (n<e-^ t )e- B "V Ql ,.., Qjv (x;t), 

i=l 

with * as Bethe wave function 

*(x) =J2 A ° eMp) '*- ( 22 ) 

(7 

^ is a tensor of rank N with components ^ Ql ... ajv (x), where the sum is taken 
over all permutations a of (1, 2, 3, N). Inserting (21) in (4) therefor the above 
eigenfunctions correspond to the eigenvalues as follows 

N 

i=i 

The next step is to determine the coefficients A CT . Inserting (21) in the type 1 
boundary condition (6) gives 

(24) 

relation (24) is the boundary condition for ^ but it docs not has the form of 
boundary condition (6) (because of the time dependency and the elements of 
power x). To build a modified form of (24) (similar to (6)), we should omit 
the time dependency and the elements of power x from both-hand sides of this 
equation. After finding a modified form, we can put (22) in it to obtain the 
coefficients A CT . Obviously to obtain a modified desired form one arrives at 

v ai + v a2 = vp 1 + vp 2 , (25) 

v ai v a2 = v fjl v f 3 2 . (26) 
one can write these equations as 

v ai + v a2 = Vfa + v 02 , (27) 



(21) 



(23) 
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<+vl 2 =v} i+ v%. (28) 

These equations have a physical interpretation in terms of an elastic collision. In 
the elastic collision of equal masses (supposing that vp t s are the velocities before 
collision and v ai s are the velocities after collision) the conservation of momentum 
and kinetic energy are the same as (27) and (28) respectively. Solving these 
simultaneous equations we get 

{v ai = v^,v a2 = vp 2 } and {v ai = Vp 2 ,v a2 = v^} (29) 

Two above solutions imply that only two states for (3 are (3 = (0:1,0:2) and 
P — (02,01). These states introduce a system with PDHR and the exchange- 
reaction processes as we will show later on. So the matrix c has non-zero 
elements provided /3 = (a l7 o 2 ) and (3 — (o 2 ,Oi) are our states and therefore 
(24) reduces to 

ip ai a 2 (x,x) = cl\ll^ aia2 {x,x + l) + {^)cl\l\^ a2ai {x,x + l) 

v a 2 

= a + f3> 1 di\i\^^ 2 (x,x + l) (30) 



where 

ai<«2 

in the same way, for the type 2 boundary condition (7), one arrives at 

In the p-species case, c(c) has pl(p — 1) non-diagonal non-zero elements that 
these elements are the exchange-rates of the p-species particles but 6(6) has 
p\(p — 1) + p 2 non-zero elements (p 2 diagonal and p\(p — 1) non-diagonal) that 
describe the rates of the drop-push interactions. One should note that the 
annihilation rates are not the original elements of the matrix c(6) and they 
enter the matrix when we impose the constraint (17) (constraint (18)) to this 
matrix. Now we can write (24) in the compact notation as 

*(• • • , x k = x, x k+1 = x ,---) = c fe;fc+ i*(- • • ,x k =x, x k+1 = x + 1, • • • ). (33) 

where c k , k +i 

c k , k +i = 1 <8> • ■ • ® 1 ® ®1 ® • ■ • ® 1. (34) 

kM+l 

Note that however (29) leads to the modified desired form of (24) but in itself 
implies some loss of generality on the elements of the matrices b and c. Now 
the coefficients A CT can be determined by putting (22) in (33), which gives 

[1 _ e *(P*+i)c fc>fc+1 ]A„ + [1 - e^d k , k+1 ]A aak = 0. (35) 

from this one obtains 

A ^fc = S kl+l(<7(Pk),0-(Pk+l))Acr, (36) 



for [3 — a 
for f3 = (o 2 , oi) 



(31) 



I 1 \ljaia?, 

V al V a2 ' «lO!2 



for j3 — a 



01 



for 



(02,01). 



(32) 
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where the matrix is defined through 

S { - 1 \z 1 ,z 2 ) = -(1 - zic)-\\ - z 2 ~c). (37) 

and the definition of S^\ +1 is similar to that of Ck,k+i in (34). We have also 
used Zk = e tpk . The same procedure for the type 2 boundary condition (7) 
results in 

SW(z u z 2 ) = -(1 - z 2 x l)-\l - z^b). (38) 

Equation (36) allows one to compute all the A^'s in terms of Ai (which is 
set to unit). As the generators of the permutation group satisfy <7k&k+iVk = 
Cfc+iffcCfc+i, so one also needs 

^a k a k+1 <7 k = A (Tfc+1(Tfe(Tfc+1 . (39) 

This in terms of ^-matrices becomes 

Si 2 (z 2 , z 3 )S 23 (zi, z 3 )Si 2 (zi, z 2 ) = S 23 (zx, z 2 )Si 2 (zi, z 3 )S 23 (z 2 , z 3 ). (40) 

Writing the S'-matrix as the product of the permutation matrix II and an R- 
matrix, 

Sfc,k+i =: nfe.fe+iiJfe.fe+i, (41) 
Equation (40) is transformed to 

R2 3 {z 2 , z 3 )Ri 3 (zi, z 3 )R 12 (zi, z 2 ) = R 12 (z 1 , z 2 )R 13 (z 1 ,z 3 )R 23 (z 2} z 3 ). (42) 

This is the spectral Yang-Baxter equation. The Bcthe ansatz solution exists, if 
the scattering matrix satisfies (40), in other words the matrix c in (37) and b in 
(38) are acceptable, only if the resulting S- matrices satisfy (40). This is a very 
restricted condition and needed for having the solvability. The S'-matrices (37) 
and (38) are exactly the ones considered in [13] and [14] respectively. Using 
the fact that is a binomial of degree one with respect to z 2 and is a 
binomial of degree one with respect to z-f 1 = e~ lpi , it can be shown that SYB 
equation (40) for and reduces to the non-spectral matrix equations 

C12[C12,C 23 ] = [C12,C2 3 ]C23, (43) 

and 

b 23 [b 23 ,bx2} = [&23,&12]&12- (44) 

which are the same as [13] and [14] respectively. The above equations are much 
simpler than (40). So it is far simpler to seek the solutions of these equations 
than to seek those of (40). 

4.2 Type 1 model 

The master equation (4) with the boundary condition 

v a2 P aia2 {x, x) = <&\*2 P Pih 0. x + l) (45) 


and relation (17) describe consistently the following reactions 
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A Q — > 0A a with rate v a , 

A a Ap — > ApA a with rate c^, 

A a Ap — > with rate <5 a £, 

A a Aji — > A Q with rate 7 Q( g, 

A a Aa ->■ with rate r? a/3 . (46) 

A matrix c, or the above reactions, corresponds to an exactly solvable exchange- 
reaction diffusion system on a one-dimensional lattice, provided c satisfies (17) 
and c satisfies (43) and the non-diagonal elements are nonnegative. 

4.3 Type 2 model 

The boundary condition 

v a2 P aia2 (x,x) =Y J h ti\Phfr{x - l,x) (47) 
and (18) with the master equation (4), introduce the reactions as follows 



A a % — > 0A Q with rate v a , 

A a Af}% -> 0A,g,4 Q with rate 6^, 

A a Aj3% -> A/30 with rate <5 Q ,g, 

A a Af}% -> A Q with rate 7 a/3 , 

ylaA/30 -> with rate ry ai g. (48) 

In brief, one can find the solvable systems corresponding to 6 if the matrix 6 
satisfies a constraint (18) and 6 satisfies (44), and all the permissible elements 
of b (including diagonal) are nonnegative. 

For more-than-two particle, now we consider a system consisting N particles 
of various species, with the evolution equation 

P(xi, . . . ,x N ) = v ai P(xi-l,...,x N )-\ h v aN P(xi, . . . ,x N - 1) 

- (v ai + ... + v aN )P{xi + ... + . x N ) (49) 

in the whole physical region, and boundary condition 

P(...,x k = x, x k+ i =x,...) = 6fe,fe+iP(. ..,Xk=x- l,x k+ i =x,...), (50) 

where bk,k+i is defined like (34) and P is an TV-tensor that the components of 
which are probabilities. It is seen that in this system, apart from the simple 
diffusion, there is a reaction between a block of n + 1 adjacent particles 

A ao ■ ■ ■ A a J — -> 0A 7O ■ • • A ln with rate (6„_i,„ • • • 6o,i)° .7.'° n " , (51) 

this comes from the fact that 

P(x_ = x, . . .jX n -i = x + n-l,x n =x + n-l) 

= (6„_i, n • ■■b . 1 )P(x = x - 1, . . .,x n -i =x + n-2,x n =x + n-l). 

(52) 

Note the order of the matrices b. This order suggests that if a collection of n + 1 
particles are adjacent, there is a probability that the first particle pushes the 
second and changes the type of second (and itself) and then it is the second 
(modified) particle that interacts with the third. 
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5 Two-particle system and exact solution 



Let us focus on the type 1 model. We consider p — 2 case, but the argument 
can be easily applied to arbitrary p. Taking A\ — A and A 2 = B. One may 
label the two-particle states (ai,a 2 ) as following 



|1>=(1,1), |2 >= (1,2) , |3 >= (2, 1) , |4>=(2,2). 
the most general matrix c that satisfies (17), can be written as 



c = 



vi - Ai 








vi- X 2 - c 32 
C32 






C23 

V 2 - A3 - C 23 





(53) 



(54) 



and therefore c is (see (31)) 



/ vi - Ai 




Vl 



V 





A 2 

(Sr)c3 2 




c 32 



V2 



(£)<*3 
A 3 





C23 



v 2 - A 4 / 



(55) 



vi - A 4 / 
Now, we must check that under what 



where, according to (20), Ai = 6i+ji + r]., 
conditions, c satisfies the (43). If one writes (43) as RHS — LHS = 0, then 
one has 64 equations that must be solved for 16 variables. By using a symbolic 
manipulator (e.g. MAPLE), one obtains the most general solution. We express 
the class of parameters of this solution in the following two sets as the given 
parameters (Gp) and the arbitrary parameters (Ap) 

Gp : {c 32 = 0,£i = (S 2 + 72 + m) - (7i + Vi)}, 

Ap : {v 1 ,v 2 ,c 23 , (li, Vi), (£2, 72, »72), (<53,73,f73), (<5 4 ,74,?74)}. (56) 

Since all reaction rates must be nonnegative, thus for <5i > 0, one should choose 



.40 


V A 


®A 


H0 


V A 


0/3 


BA 




AB 


AA 




%A 


AA 




A$ 


AA 







AB 




%B 


AB 




A% 


AB 







BA 




%A 


BA 




S0 


BA 







BB 


6 A 


%B 


BB 


74 


B% 


BB 




00 



(57) 

The model built on the reactions (57) is integrable. This is a remarkable model 
because the rate of all the interactions is arbitrary except, AA 00, that the 
rate of this process, Si, depends on some other processes. 
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Some subcases of this model: In the above model, if we abandon the 
annihilation processes (Aj =0), we have the first subcase as 



M V A %A 
B0 ^ 0i3 
BA n AB 



(58) 



this is the generalization of the reactions studied in [27] in which C23 = v 2 — v\. 
It has been showed in [26] that in the case of the single-species particles, each 
particle has a distinct hopping rate according to its site instead of its type (a 
particle in the position Xi has a hopping rate Vi). Now we consider a single- 
species system including the annihilation processes, the boundary condition is 



v 2 P(x, x) = Oi - X)P(x, x + 1) 



with 



A = 5 + 7 + 77 



(59) 



(60) 



and the reactions are 



AD ^ 9A 

AA A 9A 

AA 4 A® 

AA 4 



(61) 

1 and r] = 



the reactions (61) make the second subcase. This subcase with vi 
has been considered in [28]. 

In the 2-species case for the type 2 model, the most general matrix b that 
satisfies (18) and b that is constructed from matrix b, can be written as 



/ vi - Ai 




and 



V 



b = 









V 









— A2 — &32 


&23 


&32 


V2 - A3 - 














"1— A2— 632 


&23 


V1V2 


^1" 


&32 


"2 — A3 — 023 




V\V2 









(62) 







V2 — A 4 
^5 



«2 - A 4 / 



(63) 



/ 

The matrix 6 satisfies (44) if we have two sets of Gp and Ap as 

Gp ■ = o,^i = —(63 + J3 + V3) - (71 + m)}, 

Ap : {vi,v 2 ,b 32 , (11, m), (^2,72,r? 2 ), (53,73,%), (*4, 74, m)}- 



(64) 



This case is different from that case we considered in the type 1 model. In 
the matrix c of the type 1 model the diagonal elements do not describe any 
interactions but in this model all elements of matrix b are reaction rates. For 
example en is not a reaction rate but 611 = v\ — Ai is the reaction rate of 
AA® — > %AA. So one must choose the arbitrary parameters in such a manner 
that all elements of b (reaction rates) are nonnegative. 



12 



6 The conditional probability 



Assuming that the solvability condition (43) or (44) is satisfied. If at t = we 
have the initial configuration (yi, • ■ • , Vn), we require 

P(2;t = 0) = 5 3tif (65) 

and the probability P(x; t) becomes the conditional probability (the propagator) 
U(x;t\y;0) and thus a complete solution of the problem. It seems that the 
propagator is 

N r d N 

E/(x;i|y;0) = J[ / ff^W^ 6 "^* E A ^ (p) ' x 

where the integration region for each pi is [0,27r] and A\ = 1. In the above 
expansion /(p) is the coefficient of expansion which in the second equality we 
choose it rijLi e~ 1 ^^ , and by setting pj — > + ie, where one should 
consider the limit e — > + , (66) reproduces the required initial condition (65). 
Using this propagator, one can of course write the probability at the time t in 
terms of the initial value of probability 

|P(x;t))=£V(x;i|y;0)|P(y;0)>. (67) 

v 

Note that although and are similar to ones considered in [13, 14], but 
the propagator f/^and are different since the energy spectrum and the 
definition of the propagator of our models differ from those considered there. 
For the two-particle sector, like reactions (57), there is only one matrix in the 
expression of f/Ws (c in and b in U^). So it can be treated as a c-number. 
Following [13, 14] and [17], so by using the (66), (36) and (37), the propagator 
for the two-particle system of the type 1 model is 



£/W(x;i|y;0) = \\ r] " < ' ' /' 

3=1 J 



d V E 2 t-i(p.y) 



i(pix 1 +p 2 x 2 ) _ 



(2tt)2 

1 -e l ^c 



1 - e*Pic 



-e 



(68) 



where E 2 is obtained from (23). Using the variable z\- = e tPk , a simple contour 
integration yields 



UW{x;t\y;0) = \\r :l , 
3=1 



(xi - yi)\ (x 2 - 2/2)! 



fX2—Vl+l fXl—V2 



y . • 

jr[ (x 2 - 2/1 + ! (^i - 2/2)! V ^i - 2/2 + 1 



<? -1 + 



(69) 
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and in the same way, for two-particle system of the type 2 model one arrives at 



[/( 2 )(x;t|y;0) = \[ r[ « , 



y 

{x2 - yi)l (xi -yi- V)\ 



(xi - yi)! (x 2 - 2/2)! 
t x ^~ l r , ( . {X2-Vi)b 



-1 + 



t 



(70) 



Now we concentrate in the type 1 model or As we said, the most general 

c for the two-species system is the matrix (55) with elements that determine by 
(56). So we have 



c = 



( vt - A 2 

v\ — A 2 


\ 





(t) C 23 
V2 - A3 - C23 

w 2 - A 4 ) 



(71) 



Note that only the eigenvalues of c with modulus 1 cause singularities in A CT at 
Pj = and therefore in the integrand (66). Corresponding to matrix (71), it is 
obvious that the eigenvalues of this matrix are 



vi - A 2 , 
vi - A 2 , 
w 2 - A 3 - c 23 , 
v 2 - A 4 . 



(72) 



Since all parameters that construct eigenvalues are arbitrary, so the modulus of 
eigenvalues can be 1 or different from 1. To investigate the large-time behavior 
of the propagator it is useful to decompose the vector space on which 

c acts, into two subspaces invariant under the action of c. The first subspace 
corresponding to eigenvalues with modulus 1, and the second invariant subspace 
corresponding to eigenvalues with modulus different from 1. This decomposition 
can be done by introducing two projections Q and R, satisfying 



Q + R = 1, 
QR = RQ = 0, 
[c,Q] = [c,R]=0. 



(73) 



Q projects on the first subspace and R projects on the second. Following [13], 
we multiply by Q + R = 1: 



C/ (1) (x;t|y;0) = U [1) Q + U [1) R. 



(74) 



In the term multiplied by R, one can treat c as a number with modulus different 
from 1, Thus in this term, there is no pole in and hence in A„. So the 
integrand in (66) is nonsingular at pj = 0, which have the main contributions 
at large times. Setting pj ~ in as an approximation to arrive at 



and 



1. « (-l)W 



(75) 
(76) 
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One can also approximate E as 



£(p)«E(- 1+ ^- + f )• ( 7? ) 



So, the second term of (74) for large times results in 



U {1) R = T\ v X] ~ Vi e~ v i t — / e -[(zi-?/i-t) 2 +(z2-y2-t) 2 ]/(2t) 



which is independent of c. So at large times, the second term tends to zero 
faster than t^ 1 and only the first term of (74) survives. 

If one chooses the class of parameters in such a manner that the only eigen- 
value of c with modulus 1 is 1, then has a simple behavior for t — > oo. In 
this case, cQ = Q, and one can simplify to find 



t/«(x;t|y;0) = \[ /, " ' j 



(xi - yi)\ (x 2 - 2/2)! 



+ ^ (x 2 -yi + 0!(a;i-y2)! V + *i - 2/2 + 1/ J 

t — > 00. (79) 

This is the propagator corresponding to a single-species asymmetric simple ex- 
clusion process with PDHR multiply by Q. In fact (79) is the case of the model 
in [26], times Q. 



7 Conclusion 

We have defined two generalized totally asymmetric exclusion processes, in 
which we have p-species particles that interact and diffuse where each parti- 
cle A a has its own intrinsic hopping rates v a . We have introduced two types 
of boundary conditions in terms of two p 2 x p 2 matrices c and b respectively, 
that led to the two new models, i.e. reaction-diffusion and drop-push systems 
with PDHR. At first we have assumed that these interactions preserve the to- 
tal number of particles, so we have obtained the constraint ^2 S F^ 1 ^ = vp 1 ( 
here and elsewhere F denotes either b or c). Next we have shown that when 
one violates the conservation of particles, the mentioned constraint changes into 
^2s-^alal — v Pi — 'Vifti an d the annihilation processes add to the previous re- 
actions of new models, where A ( g 1( g 2 is the sum of all annihilation processes. We 
have also shown that only for the states /3 — a and ft — (0:2,0:1), F introduces 
the new model with PDHR and has non-zero elements. These states resulted 
in that matrix F, a version of F that constructs S'-matrix and determines the 
coefficient A„, must satisfy the non-spectral equation matrix for solvability of 
our models in the sense of the Bethe ansatz. Then we have considered the 
two-particle systems and obtained the permissable reaction rates by checking 
the solution of the non-spectral matrix equation. Finally the conditional prob- 
ability of our new models has been obtained and for the two-particle sector, 
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the exact expression and specially its large-time behavior have been calculated. 
We have shown that the propagator of 2-spccics reaction-diffusion system with 
PDHR in large times is equivalent to the propagator of the single-species asym- 
metric simple exclusion process with PDHR. 
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